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ABSTRACT 

At redshifts z^bs ^ 2, quasar absorption-line constraints on space-time variations in the 
fine-structure constant, a, rely on the comparison of Mgii and Feii transition wavelengths. 
One potentially important uncertainty is the relative abundance of Mg isotopes in the ab- 
sorbers which, if different from solar, can cause spurious shifts in the measured wavelengths 
and, therefore, a. Here we explore chemical evolution models with enhanced populations 
of intermediate-mass (IM) stars which, in their asymptotic giant branch (AGB) phase, are 
thought to be the dominant factories for heavy Mg isotopes at the low metallicities typical of 
quasar absorption systems. By design, these models partially explain recent Keck/HIRES ev- 
idence for a smaller a in z^hs < 2 absorption clouds than on Earth. However, such models also 
over-produce N, violating observed abundance trends in high-Zabs damped Lyman-a systems 
(DLAs). Our results do not support the recent claim of Ashenfelter, Mathews & Olive (2004b) 
that similar models of IM-enhanced initial mass functions (IMFs) may simultaneously explain 
the HIRES varying-a data and DLA N abundances. We explore the effect of the IM-enhanced 
model on Si, Al and P abundances, finding it to be much-less pronounced than for N. We 
also show that the '■'C/'^C ratio, as measured in absorption systems, could constitute a future 
diagnostic of non-standard models of the high-redshift IMF. 
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1 INTRODUCTION 

In the past few years, evidence has emerged that the fine-structure 
constant, a = e^jfic, may have been smaller in high-redshift quasar 
(QSO) absorption systems than the value measured today on Earth 
(e.g. Webb et al. 1999; Murohv et al. 2004, hereafter MFW04). A 
possible explanation for the lower-redshift half of this result is 
that the abundances of the heavy Mg isotopes (^'Mg and ^^Mg) 
in the absorbers are much higher, relative to that of ^''Mg, than 
solar values. Recently, Ashenfelter, Mathews & Olive 1 2004a) pro- 
posed an d expand ed upon I Ashenfelter, Mathews & Olive 2004b, 
hereafter lAMO04 a chemical evolution model with an initial mass 
function (IMF) strongly enhanced at intermediate masses (IMs), 
whereby highly super-solar abundances of heavy Mg isotopes are 
produced via asymptotic giant branch (AGB) stars. In this paper 
we explore the side-effects of this model with a view to identifying 
possible observational signatures other than increased heavy Mg 
isotope abundances. 

The paper is organised as follows. The remainder of this sec- 
tion summarizes the QSO absorption-line evidence for and against 
variations in a and describes the sensitivity of that evidence to vari- 
ations in isotopic abundances. Section|5|briefly describes the nucle- 
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osynthesis of Mg isotopes in stars of different masses. In Section|3| 
we detail our chemical evolution models and compare them with 
those used by 'AMO04'. Section|31presents our main results for the 
predicted evolution of the total and isotopic abundances of various 
elements typically observed in QSO spectra. In each case we com- 
pare these with available data from QSO absorption-line and lo- 
cal stellar and interstellar medium (ISM) studies and discuss their 
dominant uncertainties. Section |5] assesses the overall validity of 
the IM-enhanced models in light of our new results and presents 
other arguments for and against such models. Section|6|gives our 
main conclusions. 



I.I Evidence for varying a from QSO absorption systems? 

The universality and constancy of the laws of nature rely on the 
space-tim e invariance of fund ame ntal constants , such as a. There- 
fore, since lMilnd J 1 9311 1 937h and lPirad iT937l) first suggested the 
time- variation of the Newton gravitational constant, a great diver- 
sity of theoretical and experimental exploration of possible space- 
time variations in fundamental constants has been pursued (e.g. see 
review injuzan 2003). 

High resolution spectroscopy of absorption systems lying 
along the lines-of-sight to background QSOs has provided particu- 
larly interesting constraints on variations in a over large spatial and 
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temporal baselines. Early work focused on the alkali-doublet (AD) 
method: since the relative wavelength separatio n between the two 
transi tions of an AD is proportional to (e.g. iBethe & SalpeteJ 
Il977h . comparison between AD separations seen in absorption sys- 
tems with those measured in the laborator y provides a simple probe 
of a variation. Sever al authors (e.g. Varshalovich &.Potekhin^994 ; 
Cowie & Songailall 995: Varshalovich, Panchuk & I vanchiklll99€ ; 
Varshalovich. P otekhin & Ivanchik 2000) applied the AD method 
to doublets of several different ionic species (e.g. C iv, Si ii, Si rv, 
Mgn and Aim). The strongest current AD constraints come from 
analysis of many Si iv absorption systems in ^ ~ 45 000 spectra: 
Aa/a = (-0.5 ± 1.3) x 10"^ iMurohvetal. 2001 jj 21 systems; 
2.0 < Zabs < 3.0) and Aa/a = (0.15 ± 0.43) x lO"' jChand et alj 
l2004bl 15 systems; 1.6<Zata <2.9)' 

Considerable recent interest h as focused on the many- 
multiplet (MM) method introduced bv lDzuba. Flambaum & WebbI 
1 1999a b) and Webb et al. ( 19^. The MM method is a general- 
ization of the AD method, constraining changes in a by utilizing 
many observed transitions from different multiplets and different 
ions associated with each QSO absorption system. It holds many 
important advantages over the AD method, including an effective 
order-of-magnitude precision gain stemming from the large differ- 
ences in sensitivity of light (e.g. Mg, Si, Al) and heavy (e.g. Fe, 
Zn, Cr) ions to varying a. At low redshift (0.5 < Zabs S 1-8), Mg 
lines, whose red transition wavelengths (/i>2700A) are relatively 
insensitive to changes in a, act as anchors against which the larger 
expected shifts in the bluer (2300 <A< 2700 A) Fe transition wave- 
lengths can be measured. At higher Zabs, transitions from Si and Al 
provide anchor lines distributed in wavelength space amongst a va- 
riety of Cr, Fe, Ni and Zn transitions which shift by large amounts 
in both positive and negative directions as a varies. This diversity at 
high-Zabs ensures greater reliability in the face of simple systematic 
effects compared with the low-Zabs Mg/Fe systems. 

The MM method, applied to Keck/HIRES QSO absorption 
spectra, has yielded very surprising results, w i th the first tenta- 
tive evidence for a varying a by Webb et al. ( 1999) becoming 
stronger with successively larger samples iMurphv et al.. .2001ai 
IWebb et alJboOlt iMurphv. Webb & Flambaudl2003h . The most 
recent Keck/HIRES constraint comes from 143 absorption systems 
over the range 0.2 < Zabs < 4.2 ( MFW04): Aa/a = (-0.57 ± 0. 1 1) x 
10"^. This result is quite internally robust: it comprises three differ- 
ent observational samples and approximately equal low- and high- 
Zabs subsamples, all of which give consistent results. Stubbornly, it 
has also proven resistant to a range of potential instrumental and 
astrophysical systematic effects (Murphv et al. 2001b, 2 0031). 

I ntriguingly, Chand et al. l2004al) (see also ISrianand et alj 
l2004h have analysed 23 Mg/Fe absorption systems in higher signal- 
to-noise ratio (S/N) spectra /rom a different telescope and spec- 
trograph, the VLT/UVES, claiming a precise, null result over 
the range 0.4 < Zabs < 2.3 : Aa/a = (-0.06 + 0.06) x 10'^ 
lOuast. Reimers & Levshakovl ^004) and Levshak ov et alj j2004 
also find null results in individual UVES absorbers. The discrep- 
ancy between the VLT/UVES and Keck/HIRES results is yet to 
be resolved. However, it is important to note that low-order distor- 
tions of the Keck/HIRES wavelength scale, as might be expected 
from simple instrumental systematic errors, produce opposite ef- 
fects on the low- and high-Zabs samples and so cannot fully ex- 
plain the HIRES-UVES difference iMurphv et al.. 2003.. 2004) . If 



Aa/a is defined as Aa/a = (a--ao)/ao, for a- and ag the values of a in 
the absorption system(s) and in the laboratory respectively. 



the HIRES result is incorrect, the nature of the contributory sys- 
tematic errors must be subtle and somewhat conspiratorial. 



1.2 Isotopic abundance evolution? 

lMurphvetaT]j2001bl) first identified the potential systematic error 
introduced into the MM method if the relative isotopic abundances 
of crucial anchor elements like Mg and Si underwent strong cos- 
mological evolution. This is because the absorption lines of the dif- 
ferent isotopes are spaced widely enough (~ 0.5 km s"') in these 
light ions to affect the measured line centroids (Fig.Q left panel). 
IMFW04I calculated Aa/a from their HIRES absorption line data 
as a function of the assumed relative isotopic abundances of Mg 
and Si. The results are presented in the right panel of Fig.Q Note 
the relative insensitivity of Aa/a in the high-Zab.s systems to the 
Si heavy isotope ratio, (^'Si-l-^°Si)/^*Si. As noted above, this is 
expected because of the greater diversity of transitions and line- 
shifts available in the high-Zabs regime. It is therefore important 
to note that evolution in the Si isotopic abundance cannot explain 
the high-Zabs HIRES results. Until recently, no information about 
the isotopic shifts in transitions from heavier ions used at high-Zabs 
(e.g. Cm, Fell, Niii, Znii) were available. However, due to their 
higher masses, the shifts are expected to be smaller with respect to 
the sensitivity of th e lines to vary i ng a. Indeed, the recent isotopic 
shift calculations of lKozlov et alj i2004b confirm this. 

However, Fig. Q emphasises the strong dependence of 
Aa/a in the low-Zabs systems on the Mg heavy isotope ratio, 
(~''Mg+^*Mg)/^''Mg. The importance of this potential systematic 
effect therefore depends entirely on the evolution of Mg isotopic ra- 
tios in QSO absorption systems. Unfortunately, no direct measure- 
ment of (-^Mg-l--*Mg)/^''Mg in QSO absorbers is currently feasible 
due to the small separation of the isotopic absorption lines. Fortu- 
nately however, Mg is one of the few elements for which stellar 
isotopic abundances can be measured t hroug h molecular absorption 
lines, in this case transi t ions of MgH. lOav & Lambert (2000) and 
lYong, Lambert & Ivans! J2003b') have shown that observed stellar 
values of (^^Mg+^*Mg)/^'*Mg generally decrease with decreasing 
[Fe/H1, as predicted in the Galactic chemical evolution models of 
iTimmes. Wooslev & WeaveJ il995h . The (normal IMF) model we 
present in Section |3| also predicts such a decrease (Fig. |2j. This 
has been used to argue that in the low-metallicity environments of 
Mg/Fe absorbers (typically, [Zn/H] — 1.0), one should expect sub- 
solar values of (^'Mg-l-^*Mg)/^'*Mg and so, if anything, one expects 
the low-Zabs values of Aa/a to be too positive iMurphv et alj200'Tbl 
2003, 2004; Chand et al. 2004a )^ 

Contrary to this trend, lYong et alJ i2003ah have found very 
high values of (^'Mg+^*Mg)/^'*Mg for some giant stars in the glob- 
ular cluster NGC 6752, which has metallicity [Fe/H] ~ -1.6. One 
star had (25Mg+2''Mg)/2'*Mg = 0.91 (cf. [(-''Mg + -*Mg)/2''Mg]o = 
0.27). Since IM AGB stars are thought to pro duce sig nificant quan- 
tities of -^Mg and ^'Mg (see Section Yong et alj j2003j) pro- 
posed that low-metallicity, IM AGB stars may have pollu ted this 
globular cluster. This prompted Ashenfelte r et alj i2004albh to pro- 
pose a chemical evolution model with a strongly enhanced popu- 
lation of IM stars as a possible explanation for the low-Zabs HIRES 
varying-or results. In the following sections we describe in detail 
the nucleosynthesis of Mg isotopes and construct a chemical evo- 
lution model, similar to that of AMO04. to investigate the various 
effects of an IM-enhanced IMF. 
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Figur e 1. Left panel: Isotopic structures fo r the relevant Mg/Si transitions from measurements lHallstadiuslll979l IPruUinger et aljfToicl IPickering et alj 
I199m) or calculations iBerengut et"ai]l2003h . Zero velocity corresponds to the isotopic structure's centre of gravity. Right panel: Sensitivity of Aa/a mea- 
sured bv lMFWO^I to variations in Mg and Si isotopic abundance variations. Absorption systems containing no measured Si lines (i.e. most low-Zabs sys- 
tems) contribute to the square points while only systems containing Si lines and no measured Mg lines (i.e. most high-Za(,s s ystems) contribute to th e 
triangular points. The large square and triangle mark the values of Aa/a obtained with terrestrial relative isotopic abundances iRosman & Tavlo3ll998ft : 
(2**Mg : -^Mg : ^*Mg)o = (79 : 10 : 11) and (^^Si : ^'Si : ^°Si)o = (92.2 : 4.7 : 3. 1). The abundance ratio of the heavy isotopes [i.e. (-^Mg : ^^Mg) and (^'Si : ^°Si)] 
was held at the terrestiial value along the abscissa. 



2 NUCLEOSYNTHESIS OF MAGNESIUM ISOTOPES 



Massive stars culminating in Type II SNe are responsible for most 
of the Mg isotopes in the present-day Galaxy. However, there is ev- 
idence that intermediate-mass stars may dominate the production 
of the neutron-rich Mg isotopes in the metal-poor regime. Accord- 
ing to standard models of stellar nucleosynthesis, the yield of the 
heavy Mg isotopes, ^^-^Mg, scales with the initial stellar metal- 
licity. Conversely, the generation of ^'^Mg from SNe II operates 
fairly independently of initial metallicity. Since massive stars alone 
are insufficient to account for the higher than expected values of 
25,26]yfg^24]y[g (jgtected in metal-poor stars (Gay & Lambert 2000; 
Yong et al. 2003b), it has been suggested that there is a supplemen- 
tal source of the neutron-rich Mg isotopes. 

Karakas & Lattanzio (2003) have shown that ^""Mg and ^*Mg 
production is substantial in metal-poor intermediate-mass stars 
(IMSs). At low metallicities, asymptotic giant branch stars are be- 
lieved to generate -^Mg and ^*Mg from a-capture onto -^Ne trig- 
gered by He-shell thermal pulsing. Temperatures at the base of the 
convective envelope in 4—6 M© stars can be high enough to burn 
^"^Mg via hot bottom burning (HBB) as well as synthesise large 
amounts of ^'Mg and ^*Mg. Based on the recent nucleosynthetic 
calculations from Karakas & Lattanzio (2003), AGB stars have 
been shown to produce sufficient quantities of ^s.zsjyjg resolve the 
discrepancy between local stellar observations and previous model 
predictions (Fenner et al. 2003). For the present study, we have 
incorporated the same grids of low- and intermediate-mass stellar 
yields that successfully reproduced the solar neighbourhood Mg 
isotopic evolution. 



3 CHEMICAL EVOLUTION MODELS 

The distribution of elements and isotopes as a function of time and 
radius was simulated for a Milky Way-like disk galaxy by numer- 
ically solving the classic set of equations (as described by Tinsley 
1980) governing gas infall, star formation, stellar evolution and nu- 
cleosynthesis. In particular, we present predictions for (i) the so- 
lar annulus (r = 8.5 kpc) because it is the solar neighbourhood 
for which we have the most comprehensive set of empirical con- 
straints, and (ii) the outer disk (r = 16.5 kpc) because it is possible 
that many QSO absorption systems correspond to the outer regions 
of spiral disks (e.g. Dessauges-Zavadsky et al. 2004). To summarise 
the details of the model: we define tr,(r, t) as the mass surface den- 
sity of species ; at time t and radius r, and assume that its rate of 
change of is given by: 

d r"'"P d}(m) 

— crj(r,t) = I if/{r,t - T„,)Yi{m,Z(r,t - T„)) dm 

d 

o-,(r, Omfaii - Xi(r, t) i/f(r, t\ (1) 

at 

where the three terms on the right-hand side correspond to the stel- 
lar ejecta, gas infall, and star formation, respectively. The star for- 
mation rate, if/, varies with the square of the gas surface density 
in our models, consistent with the empirical Schmidt (1959) law. 
y,(ra, Z(r, t - T„,)) denotes the stellar yield of i (in mass units) from 
a star of mass m and metallicity Z(r, t-r^), if>(m) is the initial mass 
function, and X, is the mass fraction of element By definition, the 
sum of Xj over all ; is unity. The total surface mass density is iden- 
tical to the integral over the infall rate. The lower and upper stellar 
mass limits, mi„„ and m^p, are 0.08 Mq and 60 Mq, respectively, 
while T„, is the main-sequence lifetime of a star of mass m. We split 
the first term into three equations that deal separately with low- and 



4 Y. Fenner, M. T. Murphy, B. K. Gibson 




2 

1 

^ 

w -2 
o 

-3 
-4 



-3 -2 -1 





[Fe/H] 








H A »Si 






-3 -2 -1 





[Fe/H] 





-2 

[Fe/H] 



. O =°Cr 
.A 53Cr 





O «Ti 

A "Ti 

□ «T1 
■ X 50Ti 












-3 -2 -1 
[Fe/H] 





O 5<Fe 

A s'Fe 

□ 58Fe 

X '"Fe 


O 




-3 -2 -1 






[Fe/H] 




Figure 2. The evolution of isotopic ratios witli metallicity [Fe/H] for the 
elements C, Mg, Si, Ti, Cr, Fe, Ni, Zn, which are used in varying ha I a 
analyses. Isotopic abundances are shown on a logarithmic scale relative to 
the dominant isotope for that element. Different curves illustrate the pre- 
dicted behaviour of the various isotopes from a standard solar radius chem- 
ical evolution model. Symbols show the corresponding solar values. In gen- 
eral, the minor isotopes accrue more slowly in the ISM with respect to the 
dominant isotope for each element. A notable exception is Ni, for which 
60,6i,62j^j^58^j (jeclines over time. The departure of Ni from the general 
rule of thumb may reflect uncertainties in its nucleosynthesis, including an 
overproduction from SNe la, that have been pointed out in the literature 
(e.g. Timmes et al. 1995; Iwamoto et al. 1999). 



intermediate-mass stars, Type la supernova (SN) progenitors, and 
massive stars. 

Further details of the numerical code employed in this study 
can be found in Fenner & Gibson (2003), Fenner et al. (2003) and 
Fenner, Prochaska & Gibson (2004). 



3.1 Infall scheme 

We assumed that the Milky Way-like disk galaxy formed dur- 
ing two main gas accretion episodes. The first occurs on a rapid 
timescale (< 0.5 Gyr) and is associated with the formation of the 
halo and thick disk, while the second episode occurs on a longer 
timescale and fuels the formation of stars in the disk. For simplic- 
ity, we have assumed no prior metal-enrichment of the gas infalling 
onto the disk, although there is some evidence from observations 



of high- velocity clouds that gas falling into the Galaxy may con- 
tain traces of heavy elements (e.g. Wakker et al. 1999; Gibson et al. 
2001; Sembach et al. 2002). Exponentially decaying infall rates 
have been adopted, such that the evolution of total surface mass 
Cimir, t) density is given by 



dt 



: A(r)e 



-l/TH(r) 



+ B{r)e 



-('-'ilcliiy)/TD(') 



(2) 



where the infall rate coefficients A{r) and B{r) are chosen in 
order to reproduce the present-day surface mass density of the 
halo/thick disk and thin disk components, which we take to be 10 
and 45Mopc"-, respectively. The adopted timescales for the in- 
fall phases are TH = 0.1Gyr and TD = 9.0Gyr at the solar radius 
rQ = 8.5kpc. Disk formation starts after an initial delay r^eiay = 
0.5 Gyr. The 'inside-out' functional form for Toir) (Romano et al. 
2000) is adopted, whereby the timescale for disk gas accretion in- 
creases linearly with radius. The Milky Way age is taken to be 
13 Gyr. 



3.2 Initial mass function (IMF) 

The shape of the stellar initial mass function (IMF) influences the 
quantity of Galactic material locked up in stars of different masses, 
which in turn determines the rate at which different elements are 
released into the ISM. The models presented in Section|4]compare 
the Kroupa, Tout & Gilmore (1993) three-component IMF with one 
enhanced in AGB stars at low metallicities. The shape of the IM- 
enhanced IMF at time f = is illustrated in Fig. |3| The height 
and width of the AGB bump was chosen to approximate that in 
AMO04. It is centre d on 5 M^ with a narrow mass range, although 
not as narrow as the IAMOM Model 1, for which the extra AGB 
population consists almost entirely of 5±1 Mq stars (compare our 
Fig. with their figure 3). Their preferred models (i.e. Models 2 
and 3) adopt a wider IM-bump and we also favour a broader peak, 
on the grounds that it is more physically realistic^. 

In our IM-enhanced models, the amplitude of the IM-bump 
decays exponentially with increasing gas metallicity. This differs 
from AMO04, who adopted time-dependent decay. Imposing a 
metallicity- rather than time-dependent IM-enhancement was de- 
signed to reflect the different physical properties of metal-poor gas 
clouds. The decreased cooling rate and magnetic field strength in 
metal-poor material is expected to influence the mass distribution 
of newly formed stars. For consistency with AMO04 Model 1, our 
scale-factor for the exponentially decaying IM-bump was chosen 
to correspond to the predicted metallicity at the solar radius at time 

= 0.2 Gyr. This corresponds to [Z/H] 1.3, or roughly l/20th 

solar metallicity. The transformation from a time- to a metallicity- 
dependent IMS burst makes it reasonable to apply the same IMF 
formalism across all Galactic radii. Because the outer disk has 
a more protracted star formation rate and builds up metals more 
slowly than in the solar radius, the effects of an IMS burst would 
be suppressed by adopting a fixed timescale. Instead, conversion to 
metallicity-dependence prolongs the temporal duration of the IMS 
burst at outer radii and amplifies the chemical signature of an IM- 
enhanced IMF. 



^ As a check, we ran a model using a narrower and taller IM-bump, con- 
serving the mass contained in the AGB -enhancement. The final results were 
similar in both cases. 
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Figure 3. The shape of the IM-enhanced IMF used in this study at time t = 
0. The IMF is composed of an AGB-biased bump superimposed on a stan- 
dard Kroupa, Tout & Gilmore (1993) 3-component power-law. The height 
and width of the AGB bump was chosen to emulate that in AMO04. The 
amphtude of the bump decays exponentially with increasing gas metallicity. 



3.3 Stellar Yields and Lifetimes 

Low and intermediate mass stars (LIMS): For stars less massive 
than 8 Mq, we incorporated the yields from the stellar evolution 
and nucleosynthesis code described in Karakas & Lattanzio (2003), 
supplemented with yields for additional isotopes from hydrogen 
to sulfur and for the iron-peak, as well as unpublished yields for 
the metallicity Z=0.0001. The wide range of heavy elements and 
isotopes incorporated in these AGB nucleosynthesis models al- 
lows us to self-consistently predict the contribution from low- and 
intermediate-mass stars (LIMS) to the elemental and isotopic abun- 
dance of many species of interest in DLA and varying a studies. 
These include C, N, Mg, Al, Si, Ni and Fe. We also ran identi- 
cal chemical evolution models using two alternative sets of LIMS 
yields: those of (i) van den Hoek & Groenewegen (1997) and (ii) 
Marigo (2001) (using mixing length parameter a = 2 for consis- 
tency with AMO04). These two alternative inputs do not predict 
yields for elements beyond O, but they provide an indication of the 
uncertainty in the CNO yields, as will be discussed in Section|^| 
Type la supernovae (SNe la): We adopted a recalculation of 
the Thielemann, Nomoto, & Yokoi (1986) W7 model by Iwamoto 
et al. (1999) to estimate the yields from SNe la. It was assumed 
that 4% of binary systems involving intermediate and low mass 
stars result in SNe la, since this fraction provides a good fit to the 
solar neighbourhood (e.g. Alibes, Labay, & Canal 2001; Fenner & 
Gibson 2003). 

Massive stars: For stars more massive than 8-10 Mq that end 
their lives in violent supernova explosions, we implemented the 
yields from Woosley & Weaver (1995). Their yields span metallic- 
ities from zero to solar and cover the mass range 1 1-40 Mq. Since 
the upper mass limit of our stellar IMF extends beyond this mass 
range, we extrapolated the 35 and 40 Mq yields for m > 40 Mq. We 
took the lower energy A' models for stars ^ 25 M© and the higher 
energy 'B' models for heavier stars. Taking note of the suggestion 
by Timmes, Woosley, & Weaver (1995) that the WW95 mass-cuts 
may have penetrated too deeply within the iron core, we have uni- 
formly halved the iron yields from these models. In the uncertain 
mass range 8-12 Mq, we make the conservative assumption that 
these stars do not synthesise new heavy elements, but expel mate- 



rial with the same heavy element abundance pattern as the gas from 
which they were bom. 

Stellar lifetimes: We adopt metallicity-dependent main- 
sequence lifetimes calculated by Schaller et al. (1992). Although 
stars lose material over the course of their evolution via stellar 
winds and planetary nebulae, this model assumes that all the mass 
loss takes place at the end of the main-sequence phase. Our predic- 
tions are not significantly affected by this simplification. 

3.4 Comparison with Ashenfelter et al. (2004) 

We now list some key differences between our model and that of 
Ashenfelter et al. (2004b) and describe the impact on the results: 

• AGB yields: We have adopted the low- and intermediate-mass 
stellar yields from Karakas & Lattanzio (2003), with additional 
yields for isotopes up to S as well as the Fe-peak, along with unpub- 
lished yields for the metallicity Z=0.0001. Thus, the AGB yields of 
species s uch as N and the Mg isotopes are internally self-consistent. 
IAM004| added the Mg isotopic yields from Karakas & Lattanzio 
(2003) onto the Marigo (2001) C, N and O yields. Although this 
means that their N yields were not drawn from the same AGB mod- 
els as the Mg yields, the predicted evolution of N is very similar in 
either case, as discussed in Section l4.6.2l 

• Type la SNe yields: AMO04 employed a metallicity- 
dependent Type la supernovae rate from Kobayashi, Tsujimoto & 
Nomoto (2000) that prohibits the formation of SNe la below a min- 
imum [Fe/H] threshold of - 1 . 1 . This metallicity-dependence was 
postulated mostly on theoretical grounds. Observationally however, 
a/Fe ratios in dwarf galaxies (e.g. Shetrone et al. 2001) and S/Zn ra- 
tios in DLAs (Pettini, Ellison, Steidel & Bowen 1999) provide evi- 
dence for SN la activity below [Fe/H] 1 . Furthermore, the trend 

of [0/Fe] with [Fe/H] seen in local stars can be satisfactorily ex- 
plained without imposing a SN la metallicity threshold (e.g. Alibes 
et al. 2001). Thus, the present study does not impose a metallicity 
threshold and the SN la rate is calculated following the method of 
Greggio & Renzini (1983) and Matteucci & Greggio (1986). While 

a [Fe/H] 1 SN la threshold would have only a minor effect on 

the solar radius results, the outer disk remains metal-poor for at 
least several Gyr and would be more sensitive to the precise SN la 
prescription, as will be discussed in Section|31 

• Standard IMF: In recalculating the model of Timmes et al. 
(1995), AMO04 adopted the Salpeter (1955) single power-law IMF. 
We adopt a three-component Kroupa et al. (1993) law, which is 
flatter at lower masses and steeper at the high end of the IMF. This 
empirically derived function leads to a lower overall effective yield 
that is more consistent with local data (Pagel 2001). Although the 
choice of IMF does influence the evolution of the chemical species 
of interest in this investigation, the impact is not large enough to af- 
fect our final conclusions. Model 1 from iAMO04 also differs from 
our IM-enhanced model through the inclusion of an exponential 0.5 
Gyr time scale before the onset of the "normal" stellar IMF com- 
ponent. Thus, at the very earliest times, most of the star-forming 
gas ends up inside stars with m ~ 5 Mq. In contrast, we retain our 
normal IMF component at all times (akin to 'aMoP's Model 2). 

• IM-enhanced IMF: While the AMO04 IM-bump decays ex- 
ponentially on a fixed timescale (ranging from 0.2 to 0.4 Gyr de- 
pending on the model), we transformed this from a time to a metal- 
licity dependence. As discussed in Section l7!2l the assumption that 
the shape of the IMF is governed primarily by the chemical com- 
position of the star forming gas cloud is presumably more physi- 
cally justified than imposing a uniform time-dependence. Convert- 
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ing from the time to the metallicity domain allows us to apply the 
same IMF prescription to systems with various star formation his- 
tories, such as the slowly evolving outer Galactic disk, which we 
investigate in this paper. 

• Dual-phase infall: AMO04 assumed that the Galactic disk 
forms during a single phase of gas infall, in order to directly com- 
pare with the results of Timmes et al. (1995). We model the Milky 
Way formation using a dual-phase infall scheme, since this has 
been shown to provide a better match than single-infall models to 
the number distribution of G- and K-dwarfs in the solar neighbour- 
hood (e.g. Chiappini et al. 1997). We stress that our final conclu- 
sions are largely insensitive to our choice of infall scheme. 



4 RESULTS 

We now present the predicted chemical evolution for the solar ra- 
dius and the outer Galactic disk in the case of an enhanced popu- 
lation of intermediate-mass stars at low metallicity. These results 
are compared with those obtained using a normal IMF and plotted 
against observations of local stars and DLAs, where applicable. In 
all the following figures, solid and dashed lines denote the evolu- 
tion at the solar (8.5 kpc) and outer radii (16.5 kpc), respectively. 
Models adopting a normal IMF are shown with thin lines, while 
thick lines correspond to models with the IM-enhanced IMF, as il- 
lustrated in Fig. 13 

4.1 Type la supernova rate and the age-metallicity 
relationship 

The evolution of the Type la supernova rate predicted by the model 
described in Section|3|is shown in the upper panel of Fig.|4] The 
incidence of Type la SNe in the solar region (solid lines) is pre- 
dicted to have risen to a peak about 7 Gyr ago and steadily declined 
thereafter. In contrast, the predicted SN la rate in the outer disk 
(dashed lines) has continued to increase up to the present-day. This 
difference reflects the more rapid and efficient conversion of gas 
into stars in regions of higher surface density. 

A consequence of an IMF that increases the number of 
intermediate-mass stars is the birth of more Type la SNe progeni- 
tors. Type la SNe are understood to be associated with binary sys- 
tems of low- and intermediate-mass stars, in which the mass lost by 
the more evolved star is accreted by its white dwarf (WD) compan- 
ion until the WD can no longer be sustained by electron degenerate 
pressure and a violent explosion ensues. Comparing the thick and 
thin solid lines in the upper panel of Fig. it can be seen that 
the IM-enhanced IMF has no significant affect on the local solar 
neighbourhood SN la rate. In the outer disk, however, the AGB- 
enhancement elevates the SN la rate by up to a factor of six between 
0.5-1 Gyr and about a factor two at 3 Gyr before converging with 
the normal IMF model after ~6 Gyr (compare the thick and thin 
dashed lines). The different impact of the IM-enhanced IMF on the 
solar and outer radii is due to a faster build-up of metals in the solar 
radius relative to the outer disk. This causes the IM-bump to de- 
cay more quickly and leads to only a minor increase in the SN la 
rate. In contrast, the slower rate of enrichment in the outer radius 
leads to a longer-lasting IM-bump, which has a significant impact 
on the SN la rate at early times. IM-enhanced IMFs at low metallic- 
ity could leave an observable trace on the cosmic SN la rate, which 
could be used to discriminate between these models (Fields et al. 
2001). 

Type la SNe are major producers of iron and are responsible 



for 1/3-2/3 of the solar Fe content (e.g. Timmes et al. 1995). Thus, 
any increased incidence of SN la leads inevitably to the production 
of more Fe, as seen in the lower panel of Fig.0where the evolution 
of [Fe/H] is plotted. The inclusion of an enhanced IM AGB popu- 
lation increases the abundance of iron in the ISM of the outer disk 
by a factor of 2-3 between 1-3 Gyr. Once again, the solar radius 
is less sensitive to the IM-enhanced IMF because the amplitude of 
the IM-bump decays on a shorter timescale. 

As well as being SN la progenitors, AGBs leave behind white 
dwarf remnants that can be indirectly detected through microlens- 
ing experiments. We found that the IM-enhanced IMF increases the 
present-day number of WD remnants by only about 5% and 10% 
for the solar and outer radii, respectively. Thus, the number of WDs 
is not expected to be a sensitive discriminant of the different IMF 
scenarios. This echos the Galactic halo results of Gibson & Mould 
(1997) and Fields, Freese, & Grafi^ (2000), who found that stellar C 
and N abundances provide far stronger constraints than WD counts 
on any enhancement in the number of LIMS at early times. 

It should be noted that although AMO04i only predict the evo- 
lution at the solar radius, if their IM-enhanced model were applied 
to the outer disk we expect that they would not predict an increase 
in either the SN la rate or [Fe/H] because their SN la prescription 
prohibits the formation of Type la SNe in metal-poor environments. 



4.2 Magnesium isotopic ratio 

Figure|5lshows the sensitivity of (^'Mg+^*Mg)/^''Mg versus [Fe/H] 
to the choice of IMF for both the solar (8.5 kpc) and outer (16.5 
kpc) radius models. Predictions are plotted against abundances ob- 
served in nearby dwarf stars by Gay & Lambert (2000; circles) and 
Yong et al. (2003b; diamonds). It has been suggested that the pro- 
duction of heavy Mg isotopes by IMSs is needed to explain the ob- 
servations in metal-poor stars (e.g. Timmes et al. 1995; Goswami 
& Prantzos 2000; Alibes et al. 2001). Indeed, Fenner et al. (2003) 
found that IMSs were responsible for most of the heavy Mg iso- 
topes in the solar neighbourhood for [Fe/H] < -1. The inclusion 
of AGB nucleosynthesis within the framework of the standard IMF 
model provides a good match to both datasets at low [Fe/H], as 
illustrated by the thin solid line in Fig.|5| At higher [Fe/H], the nor- 
mal IMF solar model matches the Gay & Lambert (2000) data but 
not those of Yong et al. (2003b). We caution that our predictions are 
best compared against the Gay & Lambert (2000) sample because 
Yong et al. (2003b) used kinematics to preferentially select halo and 
thick disk stars. Consequently, their sample contains relatively few 
thin disk members and perhaps some stars belonging to an accreted 
component (although inspection of figure 13 from Yong et al. 2003 
indicates that the fraction of accreted stars is small if one employs 
the Gratton et al. 2003 criteria to specify the accreted component). 
An uncertainty in the ^^Mg and ^''Mg yields of ±0.2 dex is illus- 
trated for the IM-enhanced solar radius model by the shaded region 
in Fig.|5| and will be discussed further in Section l4.6. II 

Although the mean (^^Mg+-'^Mg)/-''Mg for all stars with 
[Fe/H] < -1 is about 0.15, Fig. Q demonstrates that a value 5- 
9 times larger than this is required to explain the Aa/a mea- 
sure d by MFW 04, under the assumption of a solar ^^Mg:^*Mg ra- 
tio). 'aMOO? estimated that a ratio of (25Mg+2«Mg)/2^Mg = 0.62 
would remove the need for any time-variation in a. This is lower 
than our calculated value of 1.1 ± 0.3 (see Fig.0 because theirs 
is only a rough approximation and ours is based on the QSO 
absorptio n-line spectra. However, to aid comparison, we followed 
IAMO04I in imposing an IM-enhanced IMF capable of elevating 
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Figure 4. Upper panel: Type la supernovae rate as a function of time for 
the solar neiglibourliood (solid lines) and the outer disk (dashed lines) in 
the case of a normal IMF (thin lines) versus the IM-enhanced IMF (thick 
lines). Lower panel: Trend of [Fe/H] vs time. The square shows the position 
of the sun, while lines have the same meaning as in the panel above. Due 
to a faster build-up of metals in the solar radius relative to the outer disk, 
the IM-bump decays more quickly and leads to only a minor increase in the 
SN la rate and, consequently, the [Fe/H] evolution. In contrast, the slower 
rate of enrichment in the outer radius leads to a longer-lasting IM-bump, 
which has a significant impact on the SN la rate and Fe content at early 
times. 



(^^Mg-t-^''Mg)/^*Mg to -0.62 for the solar radius model (thick solid 
line). 

The behaviour of the Mg isoto pic ratios in our IM-enhanced 
solar radius model is very similar to lAMOcSI Model 1, reaching a 
maximum of (^^Mg-t-^^'Mgy^^'Mg = 0.63 at [Fe/H] = -1.65. The 
prolonged impact of the IM-enhanced IMF in the more metal-poor 
outer disk leads to a maximum (^^MgH--^Mg)/-''Mg ratio of 0.92 at 
[Fe/H] = -1.7. This is seven times higher than the corresponding 
ratio in the case of the normal IMF. If we assume that (i) the level 
of AGB stellar enhancement is a function of gas metallicity and 
declines with the build-up of metals, and (ii) that many QSO ab- 
sorption systems consist of slowly evolving objects such as dwarfs 
or outer disks of spirals, then these findings are qualitatively consis- 
tent with a scenario in which much higher neutron-rich Mg isotopic 
abundances are found in QSO absorption systems than in local stars 
of the same metallicity. 
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Figure 5. Evolution of (^'Mg+^*Mg)/^''Mg as a function of [Fe/H]. As for 
Fig|3 the solid and dashed lines show predictions for the solar radius and 
outer disk, respectively. Thin lines denote a normal IMF model and thick 
lines represent the model with an IM-enhanced IMF. Circles and diamonds 
show stellar data from Gay & Lambert (2000) and Yong et al. (2003b), 
respectively. The square shows the solar value. The shaded area indicates a 
±0.2 dex uncertainty range for the thick solid line. 



4.3 Nitrogen 

The types of AGB stars responsible for expelling enough heavy 
Mg isotopes at early times to explain an apparent variation in a are 
also thought to be the main factories of primary nitrogen. Thus, an 
enhanced AGB population should leave a strong imprint on the N 
abundances in local stars and in DLAs. The predicted variation of 
[N/q-] with [a/H] is plotted against DLA measurements in Fig.|6| 
The model results are plotted using Si as the reference element, 
while for the Centurion et al. (2003) data, the a-element is O, Si or 
S. It is reassuring that the normal IMF solar model (thin solid line) 
passes through the DLA data before rising to roughly solar value. 
However, no single homogeneous chemical evolution model will 
account for the broad spread in the DLA data and, worryingly, the 
outer disk normal IMF model (thin dashed line) overproduces N 
with respect to the data. Since the normal IMF models (thin lines) 
already tend to produce more N than is observed in DLAs, it is 
not surprising that the IM-enhanced IMF models (thick lines) over- 
produce N by more than an order of magnitude. An indication of 
the sensitivity of these results to the stellar nucleosynthesis mod- 
els is given by comparing the thick solid line with the thin dotted 
line, which corresponds to the IM-enhanced solar radius model in 
the case of van den Hoek & Groenewegen (1997) IMS yields. Sec- 
tion l4. 6. 2l describes the uncertainties in these predictions in greater 
detail. 

Inspection of Figs.|5|&|6|reveals that the empirical constraints 
imposed by DLA nitrogen abundances are strongly violated by the 
models capable of producing sufficient '^ -^Mg relative to ^''Mg to 
mimic the variation in the a obtained by MFW04.|AMO04| encoun- 
tered the same problem with their Model 1, which they sought to 
rectify in Model 2 by increasing star formation efliiciency (SFE) by 
a factor of ~2.5. FigureQshows how our [N/o-] versus [ff/H] varies 
with a factor of 2.5 increase in SFE. The agreement with the data 
is improved, but [N/a] is still too high, except at the lowest metal- 
licities. 

The curve corresponding to our increased SFE mode l in Fig.|71 
is extremely similar to the dot dash line in figure 15 from lAMOoJ 
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Figure 6. Evolution of [N/a] as a function of [ff/H]. Dashed and solid lines 
have the same meaning as in Fig|4l Model results employ Si as the reference 
ff-element. Data points with limits and error bars show DLA measurements 
by Centurion et al. (2003), where the a-element is either 0, Si, or S. The 
thin dotted line shows the IM-enhanced solar radius model in the case of 
van den Hoek & Groenewegen (1997) IMS yields rather than Karakas & 
Lattanzio (2003), to give an indication of the uncertainties associated with 
N production in the stellar models. 

We note that their Models 1 and 2 differ not just in SFE, but in 
the location, width, amplitude, and decay timescale of the IMF 
IM-bump. Nevertheless, we wish to point out a troublesome conse- 
quence of such an increase in SFE that is robust to those differences. 
Figure|8|shows the present-day gas surface density profile expected 
in the case of normal or enhanced SFE. The standard model agrees 
with the data from Dame (1993) to within a factor of two. However, 
a factor of 2.5 increase in the SFE leads to severe gas depletion, in 
conflict with the observations. While the chemical evolution model 
that applies to QSO absorption systems need not satisfy the em- 
pirical constraints from the Milky Way, we caution that increasing 
star formation efficiency in order to alleviate overproduction of N 
might not be appropriate given the gas-rich nature of many absorp- 
tion systems. 

4.4 '^C/'^C abundance ratio 

Along with any significant contribution to the abundance of N 
and the heavy Mg isotopes in metal-poor environments, low- 
and intermediate-mass stars should leave an additional observable 
chemical signature in the form of high ''C abundance relative to 
*^C. Metal-poor AGB stars with mass ~ 4 Mq are a major source 
of ''C, produced during hot bottom burning when the CN cycle 
converts '"C into '^C. Since the important factories of '"C are less 
massive 2-3 Mq stars, the '^C/'^C ratio in the ejecta of IMSs is 
strongly mass-dependent, peaking sharply between 4 to 5 Mq. Fig- 
ure |9| shows the predicted evolution of '^C/'^C as a function of 
[Fe/H]. For both the solar and outer radii models, the introduc- 
tion of an IM-enhanced IMF increases the '^C/'^C ratio 4— 6-fold 
in the metallicity range corresponding to the greatest enhancement 
of the heavy Mg isotopes. To estimate the uncertainty in these re- 
sults due to the stellar yields, we also ran identical models using 
the LIMS yields from van den Hoek & Groenewegen (1997) and 
Marigo (2001), depicted with thin dotted and dot-dashed lines, re- 
spectively. High values of '^C/'^C may conflict with the symmetry 
of Civ line profiles observed in some QSO absorbers and could 
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Figure 7. Sensitivity of [N/o-] vs [ff/H] to the star formation efficiency. The 
solid line represents the IM-enhanced solar radius model, as seen in Figl^ 
The dotted line shows the effect of increasing star formation (SF) efficiency 
by a factor of 2.5. Increased SF efficiency improves the agreement with the 
DLA data (symbols), but still overproduces N at [a/H] — 1 . 
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Figure 8. Present-day gas surface density profile. The dotted line presents 
the result for the model with star formation efficiency increased by a factor 
of 2.5 with respect to the standard model (solid line). The dot-dashed line 
corresponds to data from Dame (1993) based on observations of atomic 
and molecular hydrogen. The standard model replicates the observations to 
within a factor of two, whereas increased star formation efficiency leads to 
excess depletion of interstellar gas. 

impose a further constraint on the shape of the IMF at early times. 
This will be discussed in more detail in Section|6| 

4.5 Silicon, aluminium and phosphorus 

Since transitions from Si are important anchor lines in the 
many-multiplet method, we show the predicted behaviour of 
(2''Si+3''Si)/2**Si as a function [Fe/H] in Fig.lTOl The relative abun- 
dance of Si isotopes is not affected by an enhanced early population 
of AGB stars to the extent of N and the heavy Mg and C isotopes. 
Indeed, (^'Si+^''Si)/^**Si remains well below the solar value (open 
square) for both the local and outer disk models (thick lines). Nev- 
ertheless, the Aff/ff measured in high-Zabs QSO absorption systems 
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Figure 9. Evolution of '""C/'^C as a function of [Fe/H] . Symbols and lines 
have the same meaning as in Fig|4] The thin dotted and dot-dashed lines 
show the IM-enhanced solar radius model in the case of van den Hoek & 
Groenewegen (1997) and Marigo (2001) IMS yields, respectively. 
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Figure 11. Variation of [Al/Fe] as a function of [Fe/H]. Dashed and 
solid lines have the same meaning as in Fig|4] DLA measurements from 
Prochaska et al. are plotted with solid square symbols, while circles and 
open squares show stellar observations from Edvardsson et al. (1993) and 
Cayrel et al. (2004), respectively. 




Figure 10. Evolution of (-''Si+3°Si)/-^Si as a function of [Fe/H]. Symbols 
and lines have the same meaning as in Fig|4l 



is largely independent of the Si isotopic composition, as highlighted 
in Section lT!^ Thus, even a dramatic elevation in (^'Si+'''Si)/''*Si 
due to AGB s tars would not account for the high-redshift results 
from llVlFWoi 

Aluminium is another element whose production by AGB 
stars is expected to be important (e.g. Karakas & Lattanzio 2003). 
The synthesis of Al via the Mg-Al chain is particularly efficient 
in metal-poor 4-6 Mq stars. The evolution of [Al/Fe] with [Fe/H] 
in the ISIVI is depicted in Fie. 1111 Observations of Al in both DLAs 
(solid squares) and nearby stars (circles and open squares) covering 
a wide range of metallicities make this element a potentially use- 
ful discriminant of the different IIVIF models. Unfortunately, most 
of the DLA measurements are lower limits, owing to line satu- 
ration (Prochaska & Wolfe 2002). The normal IIVIF models (thin 
lines) pass through the lowermost DLA detections, but only the 
IM-enhanced DVIF models (thick lines) can satisfy the highest DLA 
upper limits. The power of Al to constrain the shape of the early 
IMF is diminished by the inability of either solar radius model to 
match the stellar observations at high and low metallicity. The over- 



production of Al with respect to metal-poor halo stars and the un- 
derproduction relative to thin-disk stars is a problem of standard 
Galactic chemical evolution models that has been noted by other 
authors (e.g. Timmes et al. 1995; Alibes et al. 2001). However, we 
note that local thermodynamic equilibrium (LTE) calculations un- 
derestimate the Al abundance in metal-poor stars. The open squares 
in Fig. ll II would be about 0.6 dex higher if non-LTE effects were 
included (e.g. Gehren et al. 2004). 

The production of phosphorus from IMSs leads to a ~0.5 dex 
increase in [P/Fe] when the IM-bump is added to the standard IMF, 
as shown in Fig. 1121 Stars in the 4-6 Mq mass range are the chief 
culprits for the additional ^'P, which is understood to be generated 
through neutron-capture onto Si followed by jS-decay to ^'P (J. Lat- 
tanzio, private communication). As is the case for Al, the difference 
between the two IMF models is significant, but the DLA measure- 
ments do not allow the elimination of either scenario. In the case of 
P, the DLA data is very scant and consists mostly of upper limits 
since the Pii line is often blended in the Lyci- forest (Dessauges- 
Zavadsky et al. 2004). Nevertheless, we present our predictions for 
Al and P in the event that the growing collection of DLA measure- 
ments may eventually provide a more stringent test of the models 
presented in this paper. 

4.6 Uncertainties 

4.6.1 Magnesium 

Denissenkov & Herwig (2003) and Denissenkov & Weiss (2004) 
showed that the production of ^'Mg and ^*Mg in a typical 5 Mg 
AGB star is fairly robust to the number of thermal-pulses, the HBB 
temperature, and the third dredge-up efficiency. The final enve- 
lope abundance of ^^Mg and -*Mg in their series of models agreed 
within ~ 0. 1 dex and ~ 0.2 dex, respectively. Fenner et al. (2004) es- 
timated the sensitivity of the AGB yields of ^^Mg and ^*Mg to the 
mass-loss prescription, finding them reduced by ~ 0.2 dex when 
the Vassiliadis & Wood (1993) mass-loss law was replaced with 
a steadier rate from Reimers (1975) (with = 3.5 on the AGB). 
The treatment of convection in AGB models has recently been in- 
vestigated by Ventura (2004), who found that the yields were very 
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Figure 12. Variation of [P/Fe] as a function of [Fe/H]. Dashed and solid 
lines have the same meaning as in Fig |3 DLA measurements from 
Prochaslca et al. are plotted with square symbols, while the circle is from 
Outram, Chaffee, & Carswell (1999). 

sensitive to the adoption of standard mixing length theory (MLT) 
versus the full spectrum of turbulence (FST) model. The latter case 
results in greater mass-loss rates, shorter AGB lifetimes, less third 
dredge-up, and reduced yields of heavy elements including CNO 
and the Mg isotopes. The FST model faces serious hurdles of its 
own, since it predicts a deficit of Na in AGB ejecta and a positive 
O-Na correlation: at odds with observations of globular cluster stars 
believed to be polluted by AGBs (Ventura, D' Antona, & Mazzitelli 
2004). 

A further uncertainty relates to the presence of a '^C pocket, 
which is neglected in the AGB models implemented in this study. 
One might expect that the inclusion of this additional source of 
neutrons would increase the production of ^^Mg and ^'Mg. How- 
ever the impact of a "C pocket on nucleosynthesis in metal-poor 
IMS should be marginal in comparison to the effects of HBB (M. 
Lugaro, private communication). Despite these unknowns, the abil- 
ity of our standard Galactic chemical evolution model to reproduce 
the Mg isotopic ratios in metal-poor stars by adding in the Karakas 
& Lattanzio (2003) yields, suggests that their Mg predictions are 
reasonably accurate. 

The total production of ^'•^''Mg with respect to ^''Mg from 
Type II SNe is uncertain to a similar degree due to internal un- 
certainties in the stellar models as well as sensitivity to the shape 
and upper mass limit of the IMF. Fenner et al. (2003) found that 
a Salpeter (1955) IMF leads to 50 per cent higher present-day 
^*Mg/-''Mg ratio than in the case of a Kroupa et al. (1993) function, 
because the Salpeter law gives rise to a greater fraction of massive 
stars (see their figs. 5 and 6). The shaded region in Fig.|5|indicates 
the ± 0.2 dex uncertainty in the results due to errors in both the 
AGB and SN II contribution. 



4.6.2 Nitrogen and ^^C^-C 

To estimate the uncertainty in the predicted evolution of N and the 
''C/'"C ratio, we ran identical models using two alternative sets 
of IMS yields: those of van den Hoek & Groenewegen (1997) and 
Marigo (2001). The Marigo yields for N were in close agreement 
with those of Karakas & Lattanzio (2003) while the implementa- 
tion of the van den Hoek & Groenewegen (1997) yields produced 



even more N. The dotted line in Fig.|6|displays the results from an 
IM-enhanced solar radius model using the van den Hoek & Groe- 
newegen (1997) yields in place of Karakas et al. (i.e. compare dot- 
ted with thick solid curve). It is clear that the problem of excess 
nitrogen is further exacerbated in this case. 

The dotted and dot-dashed lines in Fig.|9|show the behaviour 
of "C/'-C assuming van den Hoek & Groenewegen (1997) and 
Marigo (2001) IMS yields, respectively. The Marigo yields gen- 
erate the highest "C/'^C peak [with ('^C/'-C),„fl., = 0.073] and 
those of van den Hoek & Groenewegen (1997) the lowest [with 
('^C/'^C),„a.v = 0.042]. The region encompassed by these two ex- 
treme models is indicative of the level of uncertainty afflicting the 
models. 



5 DISCUSSION 

5.1 Summary of Results 

iMurphv et al] <200HJ) pointed out that t^aja measurements are 
sensitive to the isotopic composition of Mg in the gas-phase of 
the low-z absorbers used in varying-a studies. Standard models 
of chemical evolution predict that the abundance of the neutron- 
rich isotopes of Mg relative to ^'*Mg decreases with lookback time, 
makin g the MFW04 t^aja result even more significant. However, 
IAMO0 4 demonstrated that a sufficiently enhanced early population 
of intermediate-mass stars can rais e (^^''Mg+^ ''Mg)/^''Mg to the su- 
persolar levels needed to render the lMFWM Ag / a result null. 

A major problem with this scenario is the almost inevitable 
overproduction of N when compared against the DLA data (Fig.|5J. 
This is because the AGB stars responsible for significant ^'-^^Mg 
production are also important sources of N. There are various ways 
to reconcile an IM-enhanced IMF with the low N abundances in 
DLAs, but none are entirely satisfactory. For instance, since nucle- 
osynthesis within very metal-poor IMSs is not well understood, the 
theoretical N yields may be overestimated. However, in this case, it 
would be difficult to explain the high N abundance found in some 
of the most iron-depleted stars in our galaxy (e.g. Norris, Beers & 
Ryan 2000; Christlieb et al. 2004). Another way to mitigate the 
problem of excess [N/a] is through increased star formation effi- 
ciency, however this leads to severe gas depletion that may be in- 
consistent with the gas-rich nature of DLAs (Fig.|S). Finally, the 
history of many DLAs may be marked by periods of supernova- 
driven outflows, but it is not clear how such galactic winds could 
preferentially remove N but not the heavy Mg isotopes. Thus, the 
observed N abundance in QSO absorbers still stands as a robust test 
of a putative IM-enhanced early IMF. 

The enhanced '''C/'^C ratio predicted by these types of chem- 
ical evolution models (Fig. |5) holds promise as a future probe 
of the AGB contrib ution to QSO absorption line abundances. 
ICarlssonet^h995h has calculated the isotopic shifts in the ubiq- 
uitous Civ /H548/1550 alkali doublet, finding the '"^C line to lie 
Av = 10.3 km s"' to the blue of the '"C line. Recent calculations 
by Berengut et al. (in preparation) confirm these large shifts to 
within 1 per cent precision and show the '^C~'"C separation to be 
Av = 5.5 km s"'. Enhancements as large as '^C/'"C a 0.1 may al- 
ready be ruled out by the symmetric C iv line profiles observed in 
R > 45 000 (FWHM < 6.7 km s ') spectra of highly ionized QSO 
absor bers with simple velocity structure (e.g. Peti tiean & Aracill 
|2004) . However, placing limits at lower "C abundances may be 
unconvincing unless the laboratory wavelengths of the C iv doublet 
transitions can be measured with better precision than 0.4km s"' 
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jGriesmann & Klindlioool) . If the calculated isotopic separations 
for transitions of other C species (e.g. C i and C ii) prove to be as 
large as those for C iv, they may provide a more reliable probe of 
the '^^C/'"C ratio in gas more closely related to DLAs and, there- 
fore, the chemical evolution models presented here. 

Given that many QSO absorbers are thought to probe slowly 
evolving systems such as dwarf galaxies and the outer regions of 
galactic disks, we simulated chemical evolution at large galacto- 
centric radii, where metal-enrichment is more gradual. We found 
that the chemical signature from an enhanced IMS population at 
low metallicities was more pronounced in the outer disk, with re- 
spect to the solar neighbourhood. This follows directly from the as- 
sumed metallicity-dependence of the IMS enhancement. Thus, our 
IM-enhanced models predict that systems with the slowest build-up 
of metals will bear the strongest signature of AGB-pollution. 

The chemical abundance constraints imposed by QSO absorp- 
tion systems do not support excess numbers of AGB stars being 
formed in the types of systems from which lMFW04l derived a vary- 
ing a. Moreover, a significantly enhanced AGB population in the 
local Galaxy is ruled out by the sub-solar ratios of the neutron-rich 
Mg isotopes to ^''Mg in nearby metal-poor stars. We now sum- 
marise other arguments for and against the IMF having a non- 
standard shape in the early universe. 



5.2 Other arguments for and against an early, IM-enhanced 
IMF 

5.2.1 Observations of deuterium, carbon and nitrogen 

Motivated by the large observed scatter and possible trends with 
metallicity of local and high-Zabs D abundance determinations. 
Fields et al. (2001) and Prantzos & Ishimaru (2001) speculated that 
D could be destroyed within stars without an accompanying in- 
crease in metals, provided that: (i) the earliest populations of stars 
were strongly enhanced in IMSs, and (ii) zero-metallicity IMSs do 
not release their synthesized C and N. Prantzos & Ishimaru (2001) 
emphasize that if the second provision is not met, C and N abun- 
dances should be highly enhanced in high-Zubs systems if their D 
has been depleted through astration alone. They note that this is 
inconsistent with the DLA data. Indeed, based on the N/a trend 
in DLAs, Prochaska et al. (2002) presented a case for a top-heavy 
IMF, with the birth of IMSs actually suppressed in many QSO ab- 
sorption systems. 

The abundance of C and N in Galactic stars and DLAs pro- 
vides one of the strongest constraints on any enhancement in the 
population of LIMS in the early Galaxy (Gibson Sl Mould 1997; 
Fields, Freese, & Graff 2000). However, this is contingent upon 
how much C and N is released by extremely metal-poor IMSs. Fu- 
jimoto et al. (2000) provide evidence favouring the suppression of 
C and N yields from zero-metallicity IMSs, whereas Chieffi et al. 
(2001) reach the opposite conclusion. Abia et al. (2001) propose 
that an early IMF peaked between 3-8 Mq helps explain the pres- 
ence of metal-poor Galactic stars with enhanced C and N abun- 
dances, presuming that LIMS release C and N in large quantities. 
Thus, the stellar models do not offer a cohesive picture of the level 
of C and N enrichment from metal-poor LIMS, nor do observa- 
tions of DLAs and metal-poor stars clarify the situation. However, 
if AGBs are responsible for -^'^*Mg enhancements large enough to 
mimic a varying a in the redshift range 0.5 < Zabs ^ 1-8, they will 
tend to have had a much higher initial metallicity than primordial. 
Thus, their yields are not subject to the very large uncertainties af- 
flicting zero-metallicity models. 



5.2.2 White dwarfs 

Evidence from microlensing surveys initially led to speculation that 
white dwarf remnants contribute up to half of the mass of the Milky 
Way halo (e.g. Bennett et al. 1996; Alcock et al. 1995). Populat- 
ing the halo with such a large fraction of white dwarfs requires 
the early IMF to be heavily biased toward stars with initial mass 
in the range 1-8 Mg. Adams & Laughlin (1996) construc ted IMFs 
capable of producing a white-dwarf dominated halo, and lAMO04l 
adopted IMF parameters that were consistent with their deduced 
constraints. However, the white dwarf contribution to the dynam- 
ical mass of the Galactic halo has since been drastically revised 
downwards to less than a few percent (Gibson & Mould 1997; 
Flynn, Holopainen & Holmberg 2003; Garcia-Berro et al. 2004; 
Lee et al. 2004), obviating the need to invoke non-standard IMFs. 

5.2.3 Theoretical star formation 

Theoretical models of star formation suggest that the IMF of pri- 
mordial gas would be biased toward higher mass stars (e.g. Abel, 
Bryan & Norman 2000; Kroupa 2001; Hernandez & Ferrara 2001; 
Mackey, Bromm & Hernquist 2003; Clarke & Bromm 2003). This 
is due to diminished cooling efficiency in metal-poor gas. However, 
recent hydrodynamical simulations of the collapse and fragmen- 
tation of primordial gas have predicted bimodal primordial IMFs 
with peaks at about 100 and 1 M© (Nakamura & Umemura 2001) 
and 30 and 0.3 Mq (Omukai & Yoshii 2003). 

5.2.4 Extragalactic observations 

Empirical support exists for an early IMF skewed toward high stel- 
lar masses in extragalactic environments. For instance, to explain 
the high metal content of the intra-cluster medium various top- 
heavy IMFs have been proposed (e.g. Gibson & Matteucci 1997; 
Elbaz, Amaud & Vangioni-Flam 1995). The photometric properties 
of ellipticals have also been well-explained assuming a top-heavy 
IMF (Arimoto & Yoshii 1987). AMO04 note that their model is 
not inconsistent with the notion of a very early population of very 
massive stars that enrich the ISM to Z~ 10"''. 

In light of the uncertainties besetting theoretical simulations of 
star formation and the interpretation of extragalactic photometric 
data, careful analysis of the nucleosynthetic enrichment patterns 
of QSO absorption systems may be a vital step in deciphering the 
shape of the IMF in different astrophysical environments. 



6 CONCLUSIONS 

An early IMF heavily biased towards intermediate-mass stars is a 
potential explanation of the a-variation apparent in Keck/HIRES 
QSO spectra. During their AGB phase these stars produce large 
quantities of -^'^*Mg which later pollute QSO absorption systems 
at Ziibs :£ 1.8, causing an apparent shift in the Mg absorption lines 
and a spurious Aa/ff < 0. However, such a strong contribution 
from AGBs should be evident in the detailed abundance patterns 
obtained in a growing number of DLAs. Using chemical evolu- 
tion models of a Milky Way-like spiral galaxy at different radii, 
we have tested the idea that AGB stars were more prevalent in the 
metal-poor Universe. The predicted chemical ramifications of an 
enhanced AGB population are not supported by the available data, 
with the strongest argument against severe AGB pollution com- 
ing from the low [N/ff] values observed in DLAs. Contrary to the 
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conclusions of lAMO04l we do not find consistent model parame- 
ters which simultaneously explain the Keck/HIRES varying-o- and 
abundance data. We also contend that a variety of other arguments 
for invoking an IM-enhanced IMF are, at best, unclear. Future mea- 
surements of (or limits on) the '^C/'^C ratio from line-profile asym- 
metries and further DLA observations of other elements synthe- 
sised by AGB stars, such as Al and P, may provide additional con- 
straints on the IMF's shape at early epochs. 

The sensitivity of low-z Mg/Fe /S.ala measurements to iso- 
topic abundance variations in Mg (Fig. and the clear link to the 
nucleosynthetic history of QSO absorbers demonstrates the impor- 
tance of careful Aa/a measurements, even in the absence of real 
variations in a. If ha I a is measured to be zero in z^bs ^1-8 ab- 
sorbers via independent me ans (e.g. by co mbining Hi 21 -cm and 
molecular absorption lines: iDarlindllOOsI) then future samples of 
Mg/Fe absorbers might reveal Mg isotopic abundance variations in 
two ways: (i) through the bulk relative line shifts between Mg and 
Fe transitions studied previously and, (ii) if we assume that AGB 
pollution of the absorbers is not uniform, through an increased scat- 
ter in the line shifts. Note that no in creased scatter is observed in the 
current HIRES or UVES samples <MFW04lchand et alj2004^ . 
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